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MATHEMATICS.—A useful type of formula for the interpolation 
and representation of experimental results. L. H. ADams. 
Geophysical Laboratory. 


Mellor, in his Higher Mathematics,’ remarks, ‘‘The reader will 
perhaps have been impressed with the frequency with which ex- 
perimental results are referred to a series formula of the type: 

y = A+ Ba + Cz? + Da? + ---, 

in physical or chemical text books. For instance, I have counted 
over thirty examples in the first volume of Mendeléeff’s The 
Principles of Chemistry and more than this number in Preston’s 
Theory of Heat.” It is well known that for the representation 
of experimental results limited power series often leave much to 
be desired, and that such formulae adequately represent the data 
only in comparatively few instances or over short ranges; and 
yet as indicated by the above quotation this fact is seldom ap- 
preciated. It is the purpose of this note to point out the ad- 
vantage of functions other than power series and to advocate 
their more general employment in cases where power series ob- 
viously do not conform to the known general form of curve given 
by the measurements. 

By far the larger number of physical changes follow laws as 
yet unknown, and so, until the precise theoretical function is 
ascertained, the experimental results in question can be repre- 
sented only by some sort of empirical formula. In the choice of 


1 J. W. Mellor, Higher mathematics for students of physics and chemistry, p. 
273. 
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the most suitable equation to be fitted to a given set of experi- 
mental results‘ there is usually considerable latitude; but among 
the host of mathematical functions there will often be certain 
ones which from their peculiar characteristics or general shape 
are especially adapted to the data in hand. Nevertheless most 
chemists and physicists, who desire to fit their experimental re- 
sults to some equation, turn blindly to the familiar power series. 


When the graph representing the 
results is obviously non-linear, 
the quadratic or parabola, y = 
A + Br + Cr’, is first tried; if 
a satisfactory fit is not obtained 
another term is added, and so 
on until there is sufficiently close 
agreement between the calcu- 
lated and observed values thru- 
out the range of observation. 
Now it happens that in some 
cases a parabolic equation repre- 
sents the relation between two 
physical quantities with great 
exactness. A notable example is 
the resistance of pure platinum, 
which is a quadratic function 
of the temperature between 0° 
and 1000° to within 0.05 per cent 
or better. Such accidents, how- 
ever, are rare. More often than 
not a cubic or a fourth power 





/ 
/ 





Fig. 1 


equation is necessary especially for experimental observations of 
high precision. For instance, according to Bridgman* a power 
series of at least five terms—and probably more—would be re- 
quired to represent the resistance of mercury as a function of 
pressure, with an accuracy of 75 per cent, the degree of accuracy 
of the experimental work thruout the range of pressures investi- 


2 Proc. Am. Acad. 44: 237, (1909). 
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gated (0-6000 atm.) As a first step in the search for a suitable 
equation it is advantageous to plot the data and to consider the 
general form of the graph, especially its limiting characteristics. 
As an illustration let us consider the e.m.f. of a copper-con- 
stantan thermoelement as a function of temperature. If we plot 
e.m.f. against temperature we obtain a curve of the type shown 
in figure 1. 

The most important characteristic of this curve is its property 
of becoming more nearly linear the farther it departs from the 
origin. That is, the curve is asymptotic to a straight line (such 
as the dotted line shown in the figure), which does not pass 
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Fig. 2 


thru the origin. It is apparent that neither a parabola nor a 
cubic equation (nor any power series with a reasonable number 
of terms) can conform over a wide range to this essential con- 
dition, namely, that of steadily increasing linearity with increas- 
ing values of E ort. A power series can be made to approximate 
only to a greater or smaller portion of the whole curve, the de- 
gree of approximation being better the greater the number of 
terms included in the series and the smaller the portion of the 
curve dealt with. 

On the other hand there are a number of possible functions 
which have the desired general form. To derive one let us con- 
sider the course of the slope dE/dé of the curve given in figure 
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1. The derivative is finite at the origin and increases with t, 
rapidly at first and then more slowly, tending toward a constant 
value as a limit (see fig. 2). Such a curve referred to the dotted 


line as the axis is given by = = 1 — e~*, in general by « = 
b (1 — e“), C being a negative quantity ;* or when referred to the 
t — axis by equation a =a+b(1—e“). Integrating and ap- 
plying the initial condition that E = 0, t = 0, we have, E = At+ B 
(1-—e*). Where A =a+),B = a This equation is of the 


desired form, and, as we have found by trial, is highly satisfactory 
in a number of cases where the graph has the general form of 
figure I, as for instance, the relations between temperature and 
e.m.f. of copper-constantan thermoelements and between pres- 
sure and resistance of manganin or of ‘‘therlo”’ wire. 

The calculation of the constants of this type of equation offers 
no especial difficulty: in fact it involves no more time or labor 
than in the case of the cubic equation: 

E=A%t+B°+C# 

a form of equation which is often used and is in general of less 
utility. The method of evaluating the three constants, A, B 
and C is as follows: Write down three simultaneous equations con- 
taining pairs of corresponding values of E and ¢. Eliminate A 
by addition or subtraction and B by division; we then obtain 
an equation in C of the form: 

(1-—e°)—P(1-e") E,-PE, -K 

(1 —e°)-Q(l-—e“) E,-QE, 
where P = i and Q = 2. Assume now several values of C and 


with the aid of tables of the exponential compute the corre- 
sponding values of K. Then plot K against C and determine the 
proper value of C by graphical interpolation. The value of C 
being now known, B and finally A are readily calculated. 


dE 
3 Another suitable function is = + b tanh ct which would lead to the equa- 


tion FE = aT + b log cosh ct. 
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The result of the application of this equation to the copper- 
constantan element is shown in table I. The constantan wire 
consisted of selected lengths cut from a spool of No. 30 B. and 8. 
“Tdeal” wire made by the Electrical Alloy Co. The points at 
100°, 217°95 and 305°9 were obtained by immersing the ele- 
ment (enclosed in a glass tube) in the vapors of boiling water, 
naphthalene and benzophenone;‘ those at the lower temperatures 
by comparison with an element® previously calibrated against a 
resistance thermometer. As table I shows, the curve when 
passed thru the three upper points fits the three lower points 
very satisfactorily; indeed the divergence is no greater than the 
experimental error. 

TABLE I 


CoMPARISON OF OBSERVED E. M. F. (E) oF CoOPPER-CONSTANTAN ELEMENT WITH 
THAT CALCULATED FROM THE EQUATION E = 74.672 t — 13892 (1—e~ 9.00261 ¢) 


DIFFERENCE! 


TEMP. E (oBs.) MICROVOLTS E (cALc.) MICROVOLTS | MIOROVOLES 
24°30 961.1 960.9 0.2 
49°59 2010.9 | 2010.8 0.1 
78°75 3291 .7 3291 .4 0.3 
100°00 4276. (4276 .) 

217°95 10248 . (10248 .) 


305°9 15203 . (152083 .) 


11 microvolt corresponds to about 0°02. 


The above equation as compared with a series formula con- 
taining a like number of constants (i.e., a cubic equation) pos- 
sesses the following advantages, (1) In accord with what might 
be expected from the general shape of the curve a better agree- 
ment has been obtained between “observed”’ and ‘‘calculated”’ 
values. For instance if a cubic should be substituted for the 
exponential equation described, the differences in the fourth 
column of table I would each be greater than 5 microvolts 
(O°1). (2) Since values of (1 — e-*) may be obtained by inspec- 

4 Proper correction being made for barometric pressure. Cf. Adams & Johns- 
ton, Am. J. Sci., 33: 538 et seq. 1912. 

5 This element belongs to Dr. W. P. White of this Laboratory. For descrip- 


tion of the comparison of this element with the resistance thermometer, see Phys. 
Rev. 31, 159. 1910. 











474 PALACHE AND SCHALLER: HODGKINSONITE 


tion and in one operation from tables of the Descending Expo- 
nential,* the operation of calculating a series of values from the 
exponential equation involves less labor than from a cubic equa- 
tion. (3) Altho extrapolation is, in general, an operation beset 
with many pitfalls, nevertheless, if an extrapolation is unavoid- 
able, it is more reasonable to employ a function of proper general 
form than to take some random function which can not be ex- 
pected to fit the experimental results except over a small part 
of the range of observation. 


MINERALOGY.—Hodgkinsonite, a new mineral from Franklin 
Furnace, N.J. C. Pauacue, Harvard University, and W. T. 
ScHa.uer, Geological Survey. 

‘The mineral here described was sent to the Harvard Miner- 
alogical Museum for identification in April of this year by Mr. 
J. J. McGovern of Franklin Furnace, New Jersey. On being 
informed that the mineral was of a new species, material for 
further study and for analysis was freely supplied by Mr. H. H. 
Hodgkinson, M. E., Assistant Underground Superintendent of 
the mine, who first found the mineral in the mine workings and 
whose name it bears. Mr. Hodgkinson states that the new min- 
eral was found in the northern part of the ore body, in that part 
of the Parker Mine formerly known as the Hamburg Mine and 
quite near the hanging wall of the west leg of the ore body, be- 
tween the 850- and 900-foot levels. The locality was marked 
by a number of slips and faults, along some of which the mineral 
occurs. It has been found in a number of specimens during the 
year but nowhere in abundance. 

Hodgkinsonite is a hydrous silicate of zinc and manganese 
crystallizing in the monoclinic system. It occurs in seams in 
massive granular ore of the typical willemite-franklinite mixture; 
the seams are generally very thin with but a film of the mineral 
which is always associated with white barite and not uncommonly 
with plates of native copper. Locally the film thickens to a 
narrow vein and then the new mineral may show individuals up 
to 2 cm. across, sharply angular in form and apparently with 


* Such, for example, as those in Becker and Van Ostrand’s hyperbolic functions 
(Smithsonian Mathematical Tables, Washington. 1909). 
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crystal faces but in reality determined in their outline almost 
wholly by the older platy barite which encloses them. The clear 
pink color and brilliant cleavage of hodgkinsonite, together with 
the snow white barite make such specimens both striking and 
attractive in appearance. One mass of ore with a surface 20 cm. 
square is at least half covered with hodgkinsonite. In one case 
only has such a vein been found in which the angular cells formed 
by the intersecting barite plates were not wholly filled by hodg- 
kinsonite so that the latter was free to develop crystal planes. 
From this specimen three crystals of good quality were detached 
and these served to establish the axial ratio of the species. 

Other crystals were found occupying cavities in thicker veins 
free from barite. These crystals, the largest 1 cm. long, were 
much affected by solution, the faces being generally dull or 
facetted. They were accompanied by black rhombohedral crys- 
tals of pyrochroite and scalenohedral crystals of calcite, both 
later in age and encrusting hodgkinsonite. The latter is im- 
planted directly on willemite or franklinite and in one specimen 
on manganese garnet. The association and mode of occurrence 
both indicate a pneumatolytic origin for the new mineral. 

Hodgkinsonite is mon clinic with normal symmetry. The 
highly perfect cleavage, rormal to the symmetry plane, has been 
taken as the basal pinacoid. The elements were calculated from 
the angles of the forms taken as (110), (011), and (221) together 
with the inclination of the cleavage (001) to the prism zone. 
These angles follow. 


(001) ; Ke ~ a} Average of 2 readings on 2 crystals. 
(110)4 4 fs ae oa } Average of 7 readings on 3 crystals. 
coy, ‘ es a } Average of 2 readings on 1 crystal. 

conny{ "4 : =. jd } Average of 4 readings on 2 crystals. 


The axial elements thus calculated are: 


Po = 0.7254 qo =1.1114 yw = 84° 33}’. 
a:b:ec = 1.539 :1 : 1.1165 B = 84° 33}’. 
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The forms observed are c(001), m(110), (011), r(221), ¢(401), 
(305), q(552), u(322). The last four of these forms were found 
only on the etched erystals, which had been measured with con- 
siderable care before the better crystals were found. These four 
high-index forms are no doubt vicinal and due to etching of the 
crystals; they are nevertheless retained in the table of measure- 
ments and a figure showing them is given since they are charac- 
teristic of most of the specimens. 

In the following table may be found all the observations made 
on these crystals, together with the calculated angles for the 


various forms. 
TABLE 1 


ANGLE TABLE FOR HODGKINSONITE 


| 


CALCULATED MEASURED LIMITS 


Sa 

SZ 

¢ p ¢ p ¢ p hohed 

c(001)....| 90°00’ | 5°274", 90°00’ 5°273’ 5°27’— 5°28’| 2 
m(110)...| 33 08 | 90 00 | 33 10 | 90 00 | 33°08’-33°11’| | 7 
s(O11)....) 4 533) 48 154) 4 52 | 48 14 3 34-6 23 | 47 28-48 47 | 12 
r(221)....| 34 49 | 69 49 | 34 48 69 48 34 45-34 51/ 69 35-69 59 | 4 
t(401)....\-90 00 | 70 28 |—90 00 | 70 09 | 70 03-70 15 | 2 
2(305).....-90 00 | 18 52 |-90 00 | 18 13 17 36-18 35 | 3 
q(552)....| 34 29 | 73 33 | 34 30 | 72 34 | 34 01-35 10! 71 45-73 23 | 4 
5 | 55 | 2 


u(322)...|\—41 47 | 56 16 |—43 1 


55 25 42 59-43 31 


The crystals are acute pyramidal in habit as shown in figure 
1, and are dominated by the equal development of unit prism 
m and pyramid r. On the etched crystals this pyramid is re- 
placed by a group of facets whose average position corresponds 
to the symbol (552). The clinodome and prism faces are always 
smooth but often dull and poorly reflecting. The faces of other 
forms are minute and of the poorest quality. Figure 2 shows 
one of the etched crystals, doubly terminated and with the small 
faces of uncertain forms characteristic only of specimens of this 
type. 

The perfect cleavage of hodgkinsonite is parallel to the basal 
pinacoid. The density is 3.91, determined by a pycnometer on 
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a gram of small fragments which were later used for analysis. 
The hardness is a very little less than 5. 

The optical characters have been only imperfectly determined. 
The optic axial plane is parallel to (010). The mean refractive 
index is 1.73, determined by the immersion method. The color 
of the mineral varies from a bright pink to a pale reddish brown, 
the luster is vitreous, the streak white. 





Fig. 1 Fig. 2 
Fig. 1. Hodgkinsonite. Common habit of implanted crystals, showing basal 
cleavage when removed from the matrix. Forms m(110), s(011), r(221). On 
etched crystals r is replaced by ¢(552) without change of habit. 
Fig. 2. Hodgkinsonite. Doubly terminated etched crystal. Forms: c(001), 
m(110), s(011), x(305), (552), w(322). 


The mineral decrepitates when held in the blowpipe flame, 
wherein it fuses readily and quietly to a brown enamel. Heated 
in a closed tube, hodgkinsonite decrepitates strongly, splitting 
up into numerous thin cleavage scales which on further heating 
yield water and become brown in color. The mineral is readily 
soluble in acid, yielding gelatinous silica. 

The chemical analysis of hodgkinsonite was made on carefully 
selected cleavage fragments of a clear pink color and gave the 
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following results. Quarter gram samples were used as the total 
amount of available material was only somewhat over a gram. 


ANALYSIS AND Ratios or HopGkinsonitTe (spy W. T. ScHALLER) 














1 | 2 3 AVERAGE RATIOS 
teresa tie | 19.92 | 19.89 | 19.77 | 19.86 | 0.331 | 1.02 
ORR pce | 20.39 | 20.97 | | 20.68 | 0.291 
es eee | 62.93 | * | | 52.93 | 0.653 sie 
Sis ior cehatacs 0.99 | 0.88 | | 0.93 | 0.017 - 
SS SE BE | 0.04 | 0.04 | 0.001 
one | 5.77t) (5.77 | 0.321 0.99 

100.21 | 








* A duplicate determination of ZnO, of which a small amount was lost, gave 
51.38 per cent. 

+t Determined directly by fusing the mineral with sodium carbonate and col- 
lecting the water in a calcium chloride tube. A determination of the loss on 
ignition, corrected for oxidation of the manganese, gave the value 4.68 per cent. 


A doubtful trace of lead was encountered but iron and chlorine 
were absent. No water was given off by the mineral at 110°. 

The ratios yield sharply the formula 3 RO.SiO,.H.O where 
R is chiefly zinc and manganese. If the manganese, calcium 
and magnesium be arbitrarily taken together then the ratio of 
MnO + CaO + MgO to ZnO is 309 to 653 or 1 : 2.11 or nearly 
1:2. The formula may then be written MnO.2Zn0.Si0..H,O 
which may, as the water is all constitutional, be interpreted as 
Mn.(ZnOH):.Si0,. 




















PROCEEDINGS OF THE ACADEMY AND AFFILIATED 
SOCIETIES 


THE ANTHROPOLOGICAL SOCIETY OF WASHINGTON 


A special meeting of the Society was held on October 28, 1913, at the 
National Museum. 

Dr. Ates Hrpuicxa addressed the Society, his subject being The 
results of the speaker’s recent trip to Peru; with remarks on the anthro- 
pological problems of Peru, illustrated with lantern slides. In 1910 the 
speaker made a brief exploratory trip in Peru, which resulted in the ac- 
quisition of some valuable data and of important skeletal collections. 
The opportunity to extend the investigations came during the early part 
of the current year, in connection with the preparation of the anthro- 
pological exhibits for the Panama-California Exposition at San Diego; 
and as a consequence three busy months were spent on the Peruvian 
coast and in certain parts of the mountain region of Peru, in explora- 
tion of the ruined cities and ancient cemeteries. The principal objects 
of the trip were, first, the mapping out as far as possible of the anthro- 
pological distribution of the prehistoric Peruvian, more particularly the 
coast people; second, the determination of the physical type of the im- 
portant Nasca group of people, which represent one of the highest 
American cultures; third, further inquiry as to man’s antiquity on the 
west coast of South America; and fourth, the extension of the speaker’s 
researches on pre-Columbian pathology. The conclusions to which the 
speaker was formerly led were in the main corroborated. In regard to 
the mountain regions much remains to be determined in the future. 
As to the pathology of the native Peruvian before contact with whites, 
the main work can perhaps be now regarded as done, or nearly so, 
altho individual variation in different morbid processes seems inexhaust- 
ible, and much in this line remains to be secured by future exploration. 
The ground. covered was extensive and the skeletal material examined 
was enormous, the selections alone filling over thirty boxes. No exca- 
vation was practiced, attention being restricted, on the coast, to the 
bones covering the surface of ancient cemeteries, exploited by the peons, 
and to burial caves and houses in the mountains. 

Since the speaker’s trip to Peru three years ago, a change for the 
worse was observed in the state of preservation of the ancient remains. 
Also, where formerly there were seemingly inexhaustible quantities of 
skeletal material there is now a dearth of it. No such collection as 
that made in 1910, when the speaker gathered 3400 important crania, 
will ever again be possible from these regions. The major part of the 
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old population of the coast region belongs to the brachycephalic type 
intimately related to the Maya Zapotec type in the north. Wherever 
they lived, these people of the Peruvian coast were wont to practice, 
more or less, the antero-posterior head deformation. Everywhere along 
the coast there are evidences of more or less admixture with a more 
oblong headed element closely related to the Aztec and Algonquin 
types of North America. As among the North American Pueblos, 
nowhere was the aboriginal Peruvian population at any time as great 
as the relatively numerous cemeteries or ruins might lead one at first 
to suppose, for these burial grounds and ruins date from different, altho 
not far distant, periods. 

The work done, while to some extent establishing a foundation, is 
merely a fair beginning. Similar investigations and collections wait 
urgently on the anthropologist in the important districts of Piura, 
Eton, and Moquegua, on the coast; in the western sierras from the 
neighborhood and latitude of Cajamarca to those of Arequipa; and in 
the eastern highlands from Tiahuanaco to Moyobamba.- The most im- 
portant problems that await solution are (1) the derivation of the Pe- 
ruvians; (2) the time of their advent into the country; (3) the extension 
and exact physical characteristics of the Aymara and Quechua; and (4) 
the genetic relations of the Peruvian to the Argentinan and Chilean 
aborigines. Besides this there remains to be established in many places 
the correlation of culture with the physical type of the people. The 
speaker repeated what he said in a former report, that, due to the lack 
of scientific supervision of a great majority of the excavations practiced 
in Peru to the present time, the archeological collections from that 
country are made up of little more than curiosities which it is impossible 
to refer either to any definite people or period. 

DantEL Foikmar, Secretary. 

















